Alternative Contributions to the Angular Correlations Observed at RHIC 

Associated with Parity Fluctuations 



Scott Pratt 

Department of Physics and Astronomy and National Superconducting Cyclotron Laboratory, 

Michigan State University 
East Lansing, Michigan 48824 
(Dated: February 9, 2010) 

Recent measurements at RHIC of angular correlations of same-sign vs. opposite sign pairs have 
been interpreted as evidence for large-scale fluctuations of parity-odd fields. In this paper, we provide 
alternative explanations of the same phenomena based on correlations from charge and momentum 
conservation overlaid with elliptic flow. These effects are shown to produce correlations with similar 
magnitudes as those measured. Other correlations are also considered, but estimates of their size 
suggest they are inconsequential. 

I. INTRODUCTION 

The prospect of observing large fluctuations of parity-odd observables in heavy ion collisions has recently 
gained great attention due to measurements by the STAR collaboration at the RHIC (The Relativistic Heavy 
Ion Collider) pQ. It has been proposed that the fluctuations derive from the color flux tubes, which carry both 
color-electric and color-magnetic flux parallel to the beam direction for the first ~ 1/2 fm/c of the collision 
[H [3]. The magnetic field lines, which have V ■ B a ^ in QCD, can be randomly parallel or anti-parallel to 
the electric field lines. However, since each flux tube might be responsible for dozens of particles, there might 
exist the opportunity to view these fluctuations. Observation of such fluctuations would validate the existence 
of coherent color magnetic flux, which is a basic feature of a non-Abelian gauge theory. 

The experimental manifestation of fluctuations comes from the parity-odd value of E a • B a , which is randomly 
positive or negative within a given tube. Via anomalous chiral couplings, these fields couple to the electromag- 
netic field EB. In the participant region of a heavy ion collision, where there exists a strong coherent transverse 
magnetic field B due to the passing ions, one can then generate a non-zero electric field due to the coupling 
with the color fields. This electric field is then randomly parallel or anti-parallel to the electric field, and should 
last a few tenths of a fm/c. Since each flux tube generates an electric field with random signs, the averaged 
electric field fluctuates as l/V^^bcs, where A/tubes is the number of such tubes. Angular correlations might 
then scale as 1/A^ u bes- Since there are fewer tubes than particles, one might hope that this fluctuation stands 
out compared to random fluctuations. 

The size of the correlations observed by STAR is of the order of 10 -3 to 1CP 4 , and given that these typically 
have multiplicities, M, of a few hundred particles at mid-rapidity, the correlations are smaller than 1/M and 
much smaller than l/A^bes; so one cannot neglect correlations that involve a handful of particles, such as charge 
or momentum conservation. The manifestations of charge and momentum conservation are discussed here, and 
are shown to provide effects which should explain the bulk, though perhaps not all, of the observed correlation. 
Other contributions are also discussed, but are found to be significantly smaller than the observed signal. 




FIG. 1: The overlap of two heavy ions from the perspective of looking down the beam line. In the left panel, one sees 
the magnetic field lines from the two ions, which add constructively in the participant (dark grey) region. These field 
lines couple to the E a • B a values which fluctuate randomly within the individual color flux tubes, to make electric fields 
illustrated in the right-hand panel. The signs of the E fields are of the same sign within a given flux tube, which extends 
over a large range of rapidities, but when averaged over the entire participant region, largely cancel and leave an average 
field strength lessened by a factor of 1/ V-^tubcs- 
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FIG. 2: Angular correlations as measured by STAR for both Au+Au and Cu+Cu collisions. The shaded areas reflect 
systematic uncertainties in the analysis related to the elliptic flow. 



The parity-odd observable measured by STAR is the quantity, 

J a/ 3 = (C0s(f/> Q + <j)p)) (1) 

Eigq,je^ c °s(^ + ^) 
M a M p 

where a and /3 represent either the positive or negative charges, and M a and Mp are the multiplicities of each 
type per event. If the two types are identical, the i = j term is omitted from the sum. The sum is calculated in 
each event, then averaged over many events. The angles 4>i are measured relative to the reaction plane of the 
event, as pictured in Fig. [I] The observable 7 Q /3 might have been defined as simply (sin</> a sini^g), rather than 
with (cos(</> Q + 4>p)) = (cos <p a cos (f>(3 — sin0 Q sinc^a). By taking the difference with the cos (j) a cos (ftp terms, it 
would be expected that certain contributions, such as those from resonances, might be subtracted away. 

STAR's measurements for both same-sign and opposite-sign correlations are shown in Fig. [2] for Au+Au and 
Cu+Cu collisions with beams of 10CL4 GeV. If the source of the correlation was solely the induced electric field, 

the opposite-sign correlation, 7-1 would have been positive and the same-sign correlation, 7 SS would have been 

equal and opposite. Indeed, "f ss is lower than 7_| , though they are not equal and opposite. The correlations 

fall with multiplicity, as one would expect for all such correlations. Since phenomena such as jets or momentum 
conservation might affect both same-sign and opposite-sign pairs, it seems reasonable to focus on the difference, 
7_i — 7 SS , as well as individual quantities. 

The next two sections focus on contributions from charge conservation and momentum conservation respec- 
tively, with ther possible contributions being discussed in the subsequent section. Quantitative comparisons are 
presented in the final section along with conclusions. 



II. CHARGE CONSERVATION 



Whenever a positive charge is created, a negative charge arises from the same point in space time. Given the 
short mean free paths of particles in heavy ion collisions, both particles then tend to be focused in the same 
rapidity and azimuthal angle by collective flow. This results in a correlation between positive and negative 
particles, i.e., for every positive particle emitted at an angle (/>+, there tends to be a negative particle emitted 
with 0_ « <j) + and with similar rapidity. Charge balance functions [4-6J represent a measure of such correlations, 
and have already been measured as a function of relative rapidity for identified particles and for relative pseudo- 
rapidity 77 for non-identified particles. Charge balance functions are calculated by taking like-sign subtractions 
of positive and negative particles, 



(2) 
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where N a b is the number of pairs per event separated by relative rapidity A77. The separation of the balancing 
charges is consistent with the charges being emitted from regions with the same collective rapidity, with the 
separation being solely driven by the final thermal motion of the decoupling source. This spread tends to 
be of the order of a half unit of rapidity. Similarly, one would expect that any charged particle would be 
accompanied by a balancing charge with a similar relative angle. The degree of correlation can be quantified 
by (cos(A$baiancc)), which would be unity if the balancing charge was always emitted exactly along the same 
direction, and would be zero if the balancing charge were to be emitted randomly. Such preliminary correlations 
have been presented by STAR [8], and show that the balancing charge tends to be emitted in within a few 
tens of degrees of its companion in central collisions. As expected, the correlations are significantly broader in 
relative angle for peripheral collisions. 

For the case of correlations to be considered here, 

Eie+je- 008(20, + (fa - &)) 

7+ - = JUMl (3) 

t £i (cos 20* (cos A^ balancc ) - sin20 i (sinA0 balancc )(4)) 
fQ Ml ' 

where it has been assumed that only one charge j balances the charge i and that its relative angle A</>b a iancc = 
4>j — <pi. Here, /q is the fraction of the charge that will be found in the sum Y]j, and can be less than unity 
due to the finite acceptance and efficiency of the detector. 

If the balancing charge were always emitted in the same direction as <pi, one could replace (cos A^baianco) = 1 
and ignore the second term in Eq. ([3]). If (cos Ac/f) ba i ance ) were assumed to be independent of the angle </>i, and 
if the balancing charge were equally like to be found at positive and negative A^baiancci Eq. ([3| simplifies to 

7+- = 2/q — ((cos A0 ba i ancc )), (4) 

where M = M + + M_, and assuming M + w M_. The double brackets emphasize that the average is also over 
the direction of 4>. More generally, the distributions of A^baiancc might depend on the direction of 4>i. Keeping 
up to second order harmonics, 

(COS A(/>balancc)(<fc) = ( (COS A0 b alanco) ) + 2u^f ) COS 2^ , (5) 
(Sin A0balance)(</>i) = 2^/ SU1 2</>i , 



V 2V = tt~ / d(j) i (smA(f) ha i ancc }((j)i)sin2(j) i 

Z7T 



1 

2^ 
1 



d(j)i (COS A^balanco) (<f>i) COS 2(j). 



»> 



With these definitions, the general form of Eq. ^ becomes 



2/q 



7 +-= M 



(V 2 ({C0S A^balanco)) + - V^J . (6) 



Since the in-plane charge-balance pairs are expected to be more tightly correlated in angle, due to more collective 

( B ) 

flow and less surface curvature for in-plane emission, should be positive. For (\>i at 45 degrees, one expects 
the balancing charge to have a slight preference to be emitted at a smaller angle due to elliptic flow which gives 
more particles for <j) — than for (j) = 90 degrees. This suggests that v\ s ' would be negative, which would 

also give a positive contribution to 7^ in Eq. Assuming that the balancing charges are highly correlated, 

i.e. ((cos A^baiancc)) ^5 1, there is little opportunity for the distribution of balancing angles to vary with (pi 

( B) ( B) FT, 

sufficiently to make v 2c or v% s very large. Assuming ((cos A^baiancc)) is fairly large, the first term in Eq. (16|) 
should dominate the other two terms. 

The fraction /q can be related to the normalization of the charge balance function, which is near 0.5 in 
STAR measurements. About half of the reduction comes from the fact that some particles within the ±1 unit 
of rapidities covered by STAR are not recorded due to finite efficiency or for having too low of a transverse 
momentum p t . The other half of the reduction comes from the balancing charge lying outside the rapidity 
range of the detector. For a particle with rapidity of zero, the balancing charge tends to be largely inside the 
acceptance, while for a particle emitted barely within the rapidity acceptance, the balancing charge has a 50% 
chance of being outside. Thus, overall the contribution from Jq due to the balancing charge being outside the 
acceptance is about 0.75. If one were to use the efficiency and acceptance corrected values of the multiplicity 
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in Eq. (|4j), rather than the number M of charges actually recorded, the fraction /q one would use should be 
on the order of 0.75. Given that collective flow accounts for about half of the final motion of the particles, one 
might expect cos A$ balancc ~ 0.5. For very peripheral collisions, one would expect a smaller value. If the charge 
balance function were measured fully as a function of 4> a and (ftp, all of the quantities on the right-hand side of 
Eq. a can be extracted from experimental analyses. 



III. MOMENTUM CONSERVATION 



Ignoring charge conservation for the moment, the implications of momentum conservation are best understood 
by considering the correlation calculated independently of the charge of the particle. In that case, 

_ J2ijtj cos + 4>j) Y^i^j ( cos 4>i cos 4>j - sin 4>i sin fc) 
7 " M(M - 1) = (Af tot )(A/ t ot-l) ' (?) 

where i and j are summed over both charged and neutral particles, i.e, from 1 to M tot = M + + Mq + M_. 
Momentum conservation implies, 

i i 

The effects of momentum conservation would be mostly model-independent if the definition of 7 was defined in 
terms of p x = p\ cos (pi instead of the angles alone. Instead, we will make a simplifying assumption that all 
particles have the same pt- This approximation should be good enough to estimate the effects of momentum 
conservation to within a few tens of percent. With this assumption, 

^cos^j = £sin<fo = 0. (9) 

i i 

Using this constraint in Eq. 

7 = -/, ^ (COS ^ 2 " Sin2 ^ ! (10) 

M tot 

where fp is the fraction of the momentum balance found within the acceptance, i.e., 

^COS0 4 = -/pcos</>,. (11) 



Again, using the definition of i>2, 



^ - f 'it- (12) 



Thus, the correlation 7 is generated by a combination of momentum conservation, which causes particles to be 
preferably generated in the opposite direction and elliptic flow, which gives more particles in the ±2 direction 
than in the ±y direction. 

Experimentally, the same-sign is defined as the average of 7++ and 7 , and assuming that the momentum 

balance is shared equally amongst the charges 

7«a = x(7+++7--) = -/i 3 TT-- ( 13 ) 

As stated above, this result was predicated on the assumption that all particles had the same pt- To remove 
this model dependence, one could consider the same quantity weighted with the transverse momenta of the 
particles, 



with the sums covering all charged particles. In that case, one comes up with the same answer as in Eq. (13 1, 
but without the model dependence, and with v' 2 being weighted with p\ , 



v. 



Ei Pit cos 2(f), 



7' = -/PT7^(l + (cos$ balancc )), v' 2 = . 2 . . (15) 



M tot Mtot {P; 
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In practice, only a subset of particles are measured. In that case some of the momentum balance comes from 
unmeasured particles and one might expect — 7 SS << V2/M. In the STAR experiment at RHIC for which these 
correlations were measured, tracks were measured for the central two units of rapidity. However, the initial 
colliding beams approached with ±5.5 units of rapidity with over half of the tracks emerging with rapidities 
outside the STAR acceptance. These particles can serve as a source of momentum, which can quench the 
momentum conservation condition, thus reduce the magnitude of 7 SS . However, the transverse momentum of 
a given track is more likely to be balanced by neighboring particles, which have similar rapidities. This is 
particularly true when considering the components of the momenta responsible for elliptic flow. For instance, 
consider a particle with final momentum transverse momentum p. Some of that momentum is generated by 
the passing nuclei, which impart a momentum k which is balanced by particles well outside the acceptance. A 
second component of the momentum q comes from interactions with neighbors with similar rapidity. For the 
track i, 

Pi = ki + qi, (16) 

with q\ = 0. Since the momenta k are generated at very early times, before particles can gain awareness of 
the elliptic anisotropy of the collision region, the vectors ki are isotropically distributed in the xy plane. When 
one calculates the correlation, 

, _Y^ijtj(^x^kx — ky^ky^) y^i^jikx <lx — ky ^qj ) Yli^tjil^ qJ — Qy\y^) 

M tot (p?) + M tot (p?) + A/ tot (p?) ~ {1) 

and given that the vectors k are isotropic and uncorrelated, 

7 = 2 n ,^2\ + ijx\ ■ ( 18 ) 



Applying momentum conservation, q\ = 0, 



M tot (pf) M tot (p t 2 ) 



Similarly, one can calculate vi with the same decomposition and find, 

, _ £*(# 2 -4 <)2 ) , JZ^fqf-kfqf) , E 4 (# 2 - Qy )2 ) 
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M tot (P 2 t) ' M tot (P 2 t) M tot (P 2 t) ( ' 

M tot (P 2 t) M tot (P 2 t) ' 



using the isotropy of k to eliminate the first term in the first line of Eq. ( 20 ) . By inspection of the last two 
expressions, one again finds 

<21) 

Thus, since elliptic flow is expected to be generated from local interactions, one expects the damping of elliptic 
flow from finite acceptance to be due to edge effects. If momentum conservation is confined to within a half unit 
of rapidity, the damping might be a few tens of percent at most for a large acceptance experiment like STAR. 

If one were to explore an observable without the cos(<pi + 4>j) weight, one would again have the negative 
contribution from momentum conservation, but it would be more strongly damped by the conservation being 
spread out over a much greater rapidity range. Furthermore, such an observable is significantly affected by 
jets or hot spots which make positive contributions to the correlations. Such observables are studied under the 
moniker of pt fluctuations. 

The contribution from charge correlation was ignored in the expression above for 7 SS , as if the two sources 
of correlation could be treated independently. However, this is not the case. For a charged particle observed 
with momentum p, one expects a particle of similar momentum but opposite charge. The same-sign correlation 
explicitly ignores the contribution of this balancing particle, but its momentum nonetheless needs to be balanced, 
and since the momentum balance was assumed to be spread equally amongst all charges, it magnifies 7 SS . The 



G 



momentum that must be balanced is increased by the factor (1 + ((cos A^balance)))- This gives the final estimate 
for the same-sign angular correlation, 

2 f v 

7ss ~ (1 + (cosA0 balance )), (22) 

where the factor of 2/3 came from assuming the charged multiplicity is two thirds the total multiplicity, M = 
2A/ tot /3. The same enhancement factor should also be applied to the expression for the p t weighted correlation, 

7ss = -^||(l + «cosA0 balance »), (23) 



J 2 



E»=i,mPm cos2^ 
M ch (p 2 t ) 



IV. OTHER POSSIBLE SOURCES OF ANGULAR CORRELATION 

The first source of correlation we consider derives from the electric field of the incoming nuclei, which can 

contribute to charge separation, i.e., 7-1 — j ss . As illustrated in Fig. [I] the magnetic field lines in the participant 

region combine constructively, whereas the electric field lines from the two nuclei would average to zero in the 
participant region. However, the location of the participants within that volume fluctuates given the random 
nature of the nucleon-nucleon collisions. Once the center-of-mass of the initial collisions strays from the center 
(x = y — 0) the averaged electric field is no longer zero. The highly Lorentz-contracted fields from the spectators 
impart an impulse Ap y to charged particles in the participant region. Due to the contraction the field decays 
within a few tenths of a fm/c, so only the charge from the initial state is affected. For the purposes of making 
a crude estimate here, we label that fraction Fq, and note that it might be at most a few tens of percent. The 
impulse such a charge i feels from a highly relativistic particle j is 

Apfa) =2a (yi ~/ J ' ) . (24) 

To estimate the effect, we considered spectators and participants to be randomly positioned in spheres of 7 
fm radius, then using a 35 mb cross section considered the net impulse imparted onto positive charges located 
according to the collision points from all the protons in the initial beams, which were assumed to move along 
straight-line trajectories. The rms Ap y for a charge in the participant region was ~ 2.5 MeV/c for central 
collisions and rose to ~ 5 MeV/c for very peripheral collisions. The induced (sin</>) for particles of a given 
charge, given the original distribution of (f> was random, is 



FoAp y 
2 (Pt)' 



(25) 



The contribution to the angular correlations is then, 



7ss - 7+ _ = (sin#«^^|. (26) 

Since Ap y is less than one percent of the mean p t , this effect comes out on the order of 10~ 5 , which is significantly 
smaller than the charge- and momentum-conservation effects described in the previous sections. However, for 
the most central collisions, where the conservation induced effects become very small, there would appear to be 
a chance that these effects could contribute. 

Although rather small, the expression above still significantly over-estimates the effect. This is because the 
initial impulse will be reduced when the initial momenta are absorbed and thermalized by the medium. Assuming 
the particles move for a time To before being absorbed, the charge separation in momentum space translates into 
a small dipole moment in the thermalized fireball. One can estimate the dipole moment by calculating how far 
the particles move before To, which might be a few tenths of a fm/c. Once that dipole moment is super-imposed 
onto a simple blast-wave picture of the final state, it again translates into a momentum anisotropy, but tends 
the value for (sine/)) tends to be lower by a factor Top t u±, m&x /(TR) than the expression above, where Uj_, ma x 
is the collective velocity at the edge of the fireball (~ 1), T is the decoupling temperature, and R is the initial 
size of the overlap region. This effectively reduces the effect by an order of magnitude. However, for higher 
Pt particles, one might expect the particles to avoid thermalization. Furthermore, for higher p t the fraction Fq 
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might be higher. Thus, it might be expected that a modest, at best, contribution from this effect might survive 
for higher p t . 

The contribution from fluctuating initial conditions should behave very similarly to the effect from parity 
fluctuations. The fluctuating initial conditions lead to electric fields which are on the order of 10% of the 
magnetic fields. It is more difficult to estimate the degree to which electric fields are generated by the anomalous 
coupling to the QCD fields and the magnetic field. Given that the coupling involves extra powers of both the 
QCD and electromagnetic coupling constants, the generated electric fields might well be smaller than those 
considered here. If such is the case, the effect of parity fluctuations would be far too small for observation. 

The most direct way to distinguish between the contributions from charge- and momentum- conservation from 
contributions from fluctuating initial conditions, or equivalently from parity fluctuations, would be to consider 
the rapidity range over which the correlations extend. For the conservation effects, the correlation should be 
confined to approximately on unit of relative rapidity, whereas the effect of fluctuation initial conditions should 
extend several units of rapidity. 

In addition to the charge separation due to the Coulomb force, one might also consider forces from the nuclear 
mean field, which due to the difference in the proton and neutron number would also be non-zero. Although 
the coupling constant for nuclear forces is larger by two orders of magnitude, the force is short range and only 
involves a few neighbors. The net result is difficult to predict, but might have a stronger impact than the electric 
force considered here. Again, the correlation might extend several units of rapidity, but the magnitude of the 
effect would scale differently with the centrality. 

A final source of correlation has to do with the two-particle correlation induced by final-state interactions. 
These correlations are driven by identical-particle interference, mutual Coulomb attraction or repulsion, or 
strong interactions between the particles. Strong interactions mainly contribute through resonant interactions, 
and thus mainly affect opposite sign correlations. However, these contributions are already included as part of the 
charge balance functions. A decaying neutral particle contributes to the charge balance function by having the 
balancing charge located in a narrow range of invariant relative momentum, which, depending on the invariant 
mass, might either broaden or narrow the charge balance function, but not change the overall normalization. 
Similarly, Coulomb interactions pull positives and negatives toward one another in momentum space, while 
pushing the positives and negatives further apart. Again, these forces do not change the normalization of the 
balance function, but can affect the width, leading to a larger value of ((cos A<^balance))- Identical particle 
statistics are mainly important for pions due to their larger phase space density. This draws same-sign pairs 
together in momentum space. Again, such correlations cannot change the fact that electric charge is locally 
conserved, so the normalization of the balance function is unaffected. Since the integrated strength of identical 
particle correlations is the average phase space- density, which averages to less than 10% over the acceptance, 
this implies that less than 10% of the strength of the balance function is pushed away from very small relative 
angle to larger relative angles, which might reduce ((cos A^baiancc)) by a few percent. 

It appears there are few effects that could modify the difference, 7_| — 7 SS . From the considerations in this 

section, it would seem that the combination 7_| — 7 SS should be dominated by charge conservation effects as 

described by Eq. Q, since competing effects would appear to be much smaller. Thus, this expression should 
be trustworthy to within the 10-20% level. It is easier to imagine effects that would contribute to 7 SS besides 
momentum conservation. The effects of jets being quenched differently between in-plane and out-of-plane is 
difficult to assess. 



V. COMPARISONS WITH DATA AND CONCLUSIONS 



Here, we replot the correlations measured by STAR in a way that better emphasizes how they might compare 

to correlations induced by momentum and charge conservation. The correlations 7_| — 7 SS and 7 SS are each 

presented in Fig. [3] after being multiplied by the multiplicity within the detector and divided by vi- The 
expressions for the correlation then become, 

^- ( 7+ - - Iss) = f Q ( ((cos A0 balancc )) + V2 ^- V2 ^), (27) 

ZV2 \ V2 V2 J 

3M (1 + ((cosA0 bal 

))) 

4^ 7ss - /P 2 ' 

The quantities on the left-hand-sides of the expressions above are plotted in Fig. [3] The multiplicity M was 
taken by assuming that the detector covered 2 units of pseudo-rapidity with the multiplicity being equal to 1.85 
times the number of participants. The factor 1.85 comes from a PHOBOS analysis that showed the multiplicity 
scaling with the number of participants to the 10% level [S]. The number of participants depends on the 
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FIG. 3: Angular correlations from STAR scaled by the multiplicity and V2- The correlations sensitive to charge 
correlation, (7+- — r y ss )M /2v2, are shown by red circles. If the angular width of the balance function is narrowly peaked, 
and if an observed particle has a high probability of having its charge-bala ncin g partner fall within the rapidity range 
of the detector, the scaled correlation should be < unity as shown in Eq. (27 1. Thus, it appears there is a possibility 

that 7-1 — "f ss is largely caused by elliptic flow combining with local charge conservation as described in Sec. [TTJ The 

scaled same-sign correlations, 3M"/ S S /4: V2, are represented by green squares. Based on the simplifying assumption that 
all particles have the same pt, Eq. (271 suggests that these correlations should be > —1 if solely driven by m omentum 
conservation. The impact of the assumption of equal p t can be tested by evaluating y' Bg described in Eq. ( 23 1. 



centrality and its extraction is modestly model dependent. Values of vi were also taken from the event-plane 
analysis of STAR [7J. 

To estimate the first term on the right-hand sides of Eq. (27 1 one must estimate the fraction of balancing 
charge falling within the rapidity range of the detector /q, and the width of the charge balance function 
((cos A0baianco))- The fraction /q of balancing charge found within the 2 units of rapidity covered by STAR can 
in principle be determined from balance function analyses. These analyses showed that the integrated strength 
of the balance function describes the chance of finding a balancing charge. This probability tended to be in 
the neighborhood of 0.5 for non-identified particles. However, some of this reduction was due to the non-ideal 
efficiency of the detector or the pt cutoffs. Since the quantity M used to scale the correlations was corrected 
for efficiency and for p t cutoffs, the fraction Jq should only refer to the fraction of balance charge that lay 
outside the rapidity acceptance of STAR, roughly ±1 units. Thus, one would expect Jq to be approximately 
0.75. Balance functions from STAR in term of the relative azimuthal angle have appeared in preliminary form 
and it appears that the widths have ((cos A</>b a iancc)) are near 0.5 within a few tens of percent. Combining with 
an estimate of Jq, it would not be surprising if the first term, /q((cos A ^balance))) in Eq. (27! is approximately 
1/3, which would explain a bit more than half of the observed correlation in Fig. [3j Since «2, s and V2, c have 
not been evaluated, even in preliminary form, one can only speculate as to whether those contributions could 
explain the remainder of the measured correlation. One does expect «2,c to be positive and V2, s to be negative, 
hence, they should help push the right-hand side toward the data, but with an unknown amount. Thus, one can 

only conclude that the majority of the correlation 7-1 — 7 SS seems to derive from combining charge conservation 

with elliptic flow, but the possibility remains that this source could fall short of the total observed correlation 
by a few tens of percent. 

If the measurement of 7_| — 7 SS turns out to be too large to be described by folding local charge conservation 

with elliptic flow, it would demand considering other sources of correlation. As shown in Sec. IV correlations 
deriving from the average electric field in the initial participant region being non-zero, tend to be too small to 
account for more than a percent of the signal. Even though the fields might give opposite impulses to positive 
and negative particles of the order of 5 MeV/c, those impulses are only given to those charges which exist in 
the first few fm/c, and then tend to be muted by the particles being reabsorbed and re-thermalized in the 
medium. The same considerations apply to electric fields generated by coupling to parity fluctuations in the 
QCD sector. If the generated electric fields are 10% of the magnetic field strength, the imparted momenta are 
also on the order of 5 MeV/c. Again, the considerations of having only a fraction of charge affected, and of 
having particles re-thermalized, would make such a signal far too small. Thus, if coupling electro-magnetic and 
QCD fields through the anomalous coupling is to account for a significant fraction of the observed correlation, 
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the generated electric field would have to be nearly as strong as the driving magnetic field. 

The correlations 7-1 and 7 SS represent moments of the full differential angular correlations which would be 



functions of <p a and (f>p. In fact, the relation for 7^ — 7 SS in Eq. (271 is more a consistency check than a test 

of charge conservation. The most telling way to distinguish correlations is to view the correlations in the most 
differential form possible, given the limited statistics. For instance, the balance function can be evaluated in 
terms of the two angles, 

S( ^ )= N + ( M + N^ a ) ' (28) 

This observable can be thought of as the conditional distribution for seeing an extra particle of opposite sign 
at fa, given the observation of a first particle at 4>\, and integrating over </>2 should yield unity for perfect 
acceptance and efficiency (also assuming equal numbers of positives and negatives). Subtractions of like-sign 
and opposite-sign correlations would carry the same information. Local charge conservation leads to a peak 
for <j>i near fa. For the parity violation hypothesis, or for the electric field acting on a fluctuating density 
distribution, one would expect the balance function to have a peak for fa — fa ~ 180 degrees, and the peak 
should be more pronounced for fa near ±90 degrees. Preliminary charge balance functions have not shown 
such a peak for A<fi ~ 180° , but analyses have not yet been performed where fa is constrained to be near ±90° . 
Furthermore, one can help discriminate sources of correlation by binning according to relative rapidity. If the 
source of the correlation is charge conservation plus elliptic flow, the peak should be confined in relative rapidity 
to one unit, the width of the measured balance function in relative rapidity. Since correlations from fluctuating 
initial conditions extend over multiple units of rapidity, one would expect the correlation due to fluctuating 
initial conditions combined with electric fields from the incoming ions, described in Sec. |IV[ to extend over a 
wide range of rapidity. Assuming that the direction of the color magnetic flux also extends over the entirety of 
a flux tube, one would also expect parity fluctuation induced correlations to also extend well beyond one unit 
of rapidity. The high statistics of recent data sets at RHIC indeed make it possible to perform such differential 
analyses. 



The right-hand side of Eq. ( 27 ) for the scaled same-sign correlation can also only be estimated at the current 
time. The fraction fp describes the fraction of the momentum responsible for V2 that falls within the same 
rapidity range must be less than unity, and if the rapidity range is dominated by final thermal motion, one 
might expect, within a few tens of percent, fp ~ 3/4. If ((cos A^baiancc)) ~ 0.5, the right-hand side would be 
approximately 1/2, which would account for only half the signal shown in Fig. [3j The estimate for j ss from 



momentum conservation in Eq. (27) was based on an a model where all the particles had the same p t . As 
emphasized in Sec. Ill evaluating the p t weighted correlation as defined in Eq. ( [23] ) should reduce the model 
dependence and give a better feci for how much of the observed correlation comes from momentum conservation 
combined with elliptic flow. It is possible that other sources of correlation are responsible for a significant 
fraction of the signal. For example, the quenching of jets, or mini-jets, might have a reaction-plane dependence. 

If the same-sign correlation is indeed largely driven by momentum conservation combined with anisotropic 
flow, this correlation might provide important insight into dynamics and perhaps even bulk properties of the 
matter. Gavin has proposed using p t correlations binned by relative rapidity as a means to gain insight into 
viscosity or heat conductivity [ID]. However, it is difficult to understand the degree to which jets interfere 
with such interpretations, as even soft-sector particles may have originated from quenched jets or mini-jets. In 
contrast, the quantity rf ss is only sensitive to the part of the momentum correlations that contribute to the 
elliptic anisotropy. Thus, if j' ss were analyzed as a function of the relative rapidity, one might gain a more 
robust insight into transport coefficients. 

Returning to the motivating suggestion that parity fluctuations drive these observables, it now seems that the 

bulk of the observed correlation 7-1 — ^ ss can be explained by the combined effects of local charge conservation 

and elliptic flow. The first term in Eq. ([6] would seem to explain half or more of the contribution, and after 
accounting for the additional terms, this has the potential to explain nearly all the signal. Further analysis of 
charge balance functions should show the degree to which the three terms in Eq. ^ contribute individually to 
the effect. The same-sign correlation, which by itself does not promote the idea of parity fluctuations, might be 
similarly explained by momentum conservation coupled to elliptic flow. However, until p t weighted correlations 
are constructed, more quantitative conclusions are unwarranted. 

Fluctuating initial conditions combined with the electric fields of the spectator nuclei in the previous section 
also provide correlations of the same sign as those observed. These effects have much in common with the effect 
one would expect from parity fluctuations. The correlations would extend multiple units of rapidity and would 
affect opposite-sign and same-sign correlations in the opposite direction. However, the estimates of the previous 
section suggest that such contributions might are too small to be observable, even though the effective electric 
fields were of the order of 10% of the strength of the magnetic fields. Given that the electric field generated 
from the magnetic field coupling to parity-odd term in the QCD fields might easily be less than 10% of the 
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magnetic field, this underscores that explaining the observed correlations through parity fluctuations would 
require surprisingly large couplings of the parity odd E ■ B terms in the QCD and electromagnetic sectors. 
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